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Abstract 

Indirect Dark Matter searches are based on the observation of secondary particles produced 

in 

, by the annihilation or decay of Dark Matter. Among them, gamma-rays are perhaps the most 

O 

<^ [ promising messengers, as they do not suffer deflection or absorption on Galactic scales, so their 

observation would directly reveal the position and the energy spectrum of the emitting source. Here, 
^ . we study the detailed gamma-ray energy spectrum of Kaluza-Klein Dark Matter in a theory with 5 

Universal Extra Dimensions. We focus in particular on the two body annihilation of Dark Matter 
particles into a photon and another particle, which produces monochromatic photons, resulting 
in a line in the energy spectrum of gamma rays. Previous calculations in the context of the five 
dimensional UED model have computed the line signal from annihilations into 77, but we extend 
these results to include 7Z and jH final states. We find that these spectral lines are subdominant 
compared to the predicted 77 signal, but they would be important as follow-up signals in the event 
of the observation of the 77 line, in order to distinguish the 5d UED model from other theoretical 
scenarios. 
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I. INTRODUCTION 



There is overwhelming evidence that our Universe contains a large component of non- 
baryonic dark matter. However, so far its identity and nature have remained elusive 
Understanding dark matter in a larger context remains one of the most compelling mys- 
teries in particle astrophysics. With the advent of the current generation of gamma ray 
and neutrino observatories, the next generation of direct detection experiments, and the 
successful operation of the Large Hadron Collider, we have entered a promising era for the 
detection of dark matter through non-gravitational interactions with the Standard Model 
(SM). 

One of the deep mysteries of dark matter is the fact that on the one hand it must be 
massive, but on the other, incredibly stable, with a lifetime on the order of the age of the 
Universe itself. This odd combination of features argues for the presence of a symmetry 
which (at least to very good approximation) forbids the dark matter from decaying. The 
most straight-forward realization is a symmetry requiring dark matter to couple in pairs 
with the Standard Model and we can classify theories of dark matter based on how they 
realize this symmetry. 

The most well known example of such a symmetry is the i?-parity often built into super- 
symmetric extensions of the SM to forbid dangerous interactions leading to large baryon- and 
lepton-number violation. While the fact that i?-parity leads naturally to a weakly interact- 
ing massive particle (WIMP) that can play the role of dark matter makes it a very attractive 
ingredient, the existence of i?-parity itself is somewhat ad hoc. It would be preferable to 
have a symmetry whose origin has a deeper motivation. 

A wide class of theories with extra dimensions naturally have such a symmetry. The 
existence of extra dimensions leads to additional spacetime symmetries, some remnant of 
which may survive compactification of the extra dimensions. In Universal Extra Dimensions 
(UED), the entire Standard Model lives equally in all dimensions. As a result, there is 
naively a translational invariance along the extra dimensional directions which translates into 

n 

conservation of Kaluza Klein (KK) number in the interactions of the modes [2J. In practice, 
the need for chiral fermions in the low energy theory selects orbifold compactifications which 
break the conservation of KK number into a Zi symmetry, but this is sufficient to guarantee 
a stable particle which can play the role of dark matter. 
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The simplest example of UED is a model with five dimensions (5d), four of which are 
ordinary spacetime with the remaining one compactified into a line segment of length L. 
The lightest Kaluza Klein (KK) mode is stable, and provided it is weakly interacting, can 
take the role of dark matter in the Universe. A standard thermal history of the Universe 

nn 

results in the correct relic density for masses on the order of several hundred GeV J3|, |4j, a 



region which is also consistent with Tevatron and precision electrowea 



discovery of at least part of the first leve 
prospects for direct |8|, |9| and indirect 10 



of KK modes at the LHC 



jam. 

£ bounds |5|, |6[ and a 
6|, [Z| • There are also 
14J detection by experiments in the near future. 



In this article, we focus in particular on the prospects for detecting gamma rays from 
annihilation of WIMPs in the galactic halo. Gamma rays offer a particularly promising av- 
enue to detect WIMP annihilations, because they point back to their sources, unlike charged 
particles, which are at the mercy of the galactic magnetic fields. However, the backgrounds 
from gamma ray sources are not very well understood, particularly around the galactic cen- 
ter where the concentration of WIMPs is greatest. Thus, we focus on the particular signal of 
a two body annihilation of two WIMPs into a photon and another particle, which produces 
monochromatic photons, resulting in a "line" in the energy spectrum of gamma rays. This 
feature is sufficiently distinct from anything produced by conventional astrophysics so as to 
perhaps balance the smaller rate caused by the fact that it is a loop-induced process 15] . 
The Fermi gamma ray observatory is actively searching for such lines as an unequivocal 



signal of dark matter [if]], and 
have appeared in the literature 17H21|. 



predictions from a wide variety of particle physics models 



Previous calculations in the context of the five dimensional UED model have computed 
the line signal from WIMPs annihilating into 77 jl^|. We extend these results to include 
7Z and final states (and we independently verify the results of pjj], finding agreement 
for the 77 channel). While both 7Z and jH turn out to be subdominant compared to 
the predicted 77 signal, they would be important as follow-up signals in the event of the 
observation of the 77 line, in order to distinguish the 5d UED model from other possibilities. 
For example, the relative strengths of the 77 and 7Z lines encodes the fact that the LKP 
couples proportionally to hypercharge, and the existence of a line, in addition to being 
interesting in its own right, is only expected to be visible for a WIMP which is a vector 
(such as the LKP) or a Dirac fermion. 

We proceed as follows. In Section [TTJ we briefly review the 5 dimensional UED model, and 
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the couplings necessary for our later computations. Sections II III and IIVI respectively cover 
the expected gamma ray continuum and lines resulting from LKP annihilations. Astrophys- 
ical inputs, the effects of finite detector energy resolution, and the gamma ray spectra are 
assembled in Section [V] We conclude in Section IVT1 



II. 5D UED MODEL 



The 5d UED model |2| consists of the ordinary four large space + time dimensions, with 
one additional spatial dimension compactified into a line segment of length L. Generic points 
in five dimensional space can be written x M = (a;**, y), where x^ is the subspace of the four 
large spacetime dimensions and < y < L is the compact dimension. The Standard Model 
field content is promoted to functions of the complete spacetime x M , with the zero modes in 
the Kaluza-Klein expansion identified with the familiar SM fields. In the discussion below, 
we will interchangeably refer to the ordinary SM fields as either "SM fields" or "zero-modes" . 
After electroweak symmetry-breaking, the lightest states obtain the usual SM masses, but 
we nonetheless continue to refer to them as zero modes in the Kaluza-Klein sense. 

Orbifold boundary conditions project out the zero modes of the fifth components of the 
gauge fields (the higher components are eaten by the KK gauge bosons as per the usual 
extra-dimensional Higgs mechanism) and the unwanted zero mode degrees of freedom of the 



fermions. The boundary conditions lead inevitably [22| to the presence of brane-localized 
terms which further shift the masses of the higher KK modes and break translational invari- 
ance along the extra dimension into a discrete parity under which the odd-numbered KK 



23 



24J . As a consequence of this Z 2 
(LKP) is stable. 
6. 



23| in which the brane 



modes are odd and the even-numbered modes are even 
remnant symmetry, the lightest first level Kaluza-Klein particle 
We follow the general picture of "minimal UED" (mUED) 
terms take sizes which are dictated by loops of bulk couplings, but we do not strictly restrict 
ourselves to the masses and couplings of mUED 1 . Consequently, the mass spectrum is chiefly 
characterized by the KK number n, with the masses of all of the particles at that level given 
roughly by 

M„ K ™ (1) 



1 For a few interesting examples of more severe deviations from mUED, see Refs. 25 1 
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The boundary terms lead to corrections for the various particles within a given level. mUED 
predicts the colored particles receive the largest positive corrections to their masses, followed 
by the SU (2)-charged states. The lightest KK modes of a given level are typically the KK 
modes of the right-handed charged leptons and the U(l) gauge boson itself. The LKP is thus 
the first level KK mode of B^, and represents a vector dark matter particle (often refered 
to as the "KK photon") which is a SM gauge singlet. Its principle interactions are to a first 
level KK fermion and its corresponding zero mode, proportional to the hypercharge of the 
fermion in question and the gauge coupling gy. 

The primary couplings of interest here are the couplings of the LKP to one KK fermion 
and one zero mode fermion, and the couplings of ordinary photons, Z bosons, and Higgs 
bosons to either zero mode or KK fermions. The photon couples to a pair of either SM or 
KK fermions proportional to eQ where Q is the electric charge as usual. The SM Higgs 
couples to a pair of zero modes proportionally to their mass, or to a pair of KK fermions 
(one of which must be an SU (2) doublet and the other an SU (2) singlet) proportionally to 
the mass of the analogue zero mode fermion. 

The SM Z boson has the standard chiral couplings to the SM fermions, and vector-like 
couplings to singlet (s) and doublet (d) KK fermions given by: 

(gv + 9a) % ( VC?) (2) 



4:S W C V 

and: 



6 (9v - 9a) %<Vtf > (3) 



As w c_ 

where £^ is the first level KK mode of a SM SU(2) doublet, and £s is the first level KK 
mode of a SM singlet. s w (c w ) are the sine (cosine) of the weak mixing angle and gy and gA 
are defined as: 

g v = -2T 3 + 4Q f s 2 w (4) 
9A = 2T 3 . (5) 

The B^ couples to one SM fermion and one KK fermion as: 

-^y^^Va + T^^ + c-c (6) 



and: 

Y "-L 
2 » 



^T^/Va-Tsj/' + cc. (7) 
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FIG. 1: The continuum gamma ray spectrum as a function of x, the ratio of the photon energy to 
the LKP mass, for LKP masses tub = 300 GeV (dash-dotted curve), 500 GeV (dashed curve), and 
1 TeV (solid curve). 

where ip^ is the zero mode fermion. 
III. CONTINUUM GAMMA RAYS 



Annihilations of two LKPs is predominantly into charged leptons (~ 59%), with signifi- 
cant fractions into quarks (~ 35%), neutrinos (~ 4%), and Higgs bosons (~ 2%) [3]]. The 
fact that a large fraction of annihilations produces hard charged leptons has the consequence 



that the continuum gamma ray signal is itself rather hard 



131 ]. which distinguishes the 



5d continuum spectrum from the six dimensional chiral square model 



17|. The resulting 



continuum spectrum is shown in Fig. [I] as a function of x = E^/Mb, the fraction of energy 
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FIG. 2: Examples of Feynman diagrams which contribute to (a) B^B^ — > 7 + Z and (b) 

of the final state photon normalized to the LKP mass. Comparing curves for LKP masses 
of 300, 500, and 1000 TeV, we see little difference in the photon spectrum, which is easily 
explained by the fact that the annihilations are almost entirely into light SM particles. 



IV. GAMMA RAY LINE CROSS SECTIONS 



In addition to a diffuse continuum of gamma rays, WIMP annihilations are also expected 
to produce prompt photons via loop-level processes. These types of annihilations produce 
7 + X final states (where X can be either a vector gauge boson or a scalar) and result 
in mono-energetic "lines" superimposed on the continuum. The energy of these lines are 
determined almost solely by the mass of the WIMP and the X particle: 

E, = M WIMP (l - Jf* ) . (8) 

Due to the non-relativistic nature of WIMPs, the possible final states (i.e., the identity of 
X) are determined by the spin of the DM particle. In the case of the 5d UED model, the 
WIMP is a vector gauge boson and, hence, X can be either a vector or a scalar. In other 
words, WIMP annihilations in the 5d UED model are capable of producing 77, 7Z and (if 
kinematically-accessible) jH final states (where H is the SM Higgs boson). 

The production of the 77 final state in the 5d UED model was first considered in Ref 
Here, we will focus on the other two possible final states (7Z and jH ) and we refer the 
interested reader to the above reference for details on the calculation of the 77 cross section. 
We have verified the results from the previous analysis and we present the numerical results 
for the flux in the following sections. 
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Representative Feynman diagrams depicting the process 

+ #(i) ^ 1 + x (where 

X = Z or H) are shown in Fig. Note that since the LKP is to very good approximation 
pure (< io^ 1 for mUED boundary terms and L^ 1 > 300 GeV 23j), loops containing W 
and/or KK charged Goldstone bosons are suppressed by the tiny content. We neglect 
these tiny contributions compared to the dominant diagrams consisting of closed loops of 
SM and KK fermions. 



A. Annihilation to 7 + Z 

For the 7 + Z process, we assign Lorentz indices and momenta as 

SffCPi) + BW(pa) 7^3) + ^ 4 (P4). (9) 



Following Ref. [IJ], the scattering amplitude involving four external gauge bosons can be 
written as: 

M = ef ( Pl )e? (p 2 )ef *(p 3 )ef *(p4)M #4lWWMl (pi,ft,J» S) P4) • (10) 

In general, the subamplitude can be expanded in terms of metric tensors and 

external momenta. Taking into account the transversality of the polarization tensors (e(p) • 
p — 0) and the non-relativistic nature of WIMPs today (such that p\ ~ p 2 = p) ; the most 
general form is given by: 

A<" 1WWW = aYaem ^ Q f \-^—P^PfP^P^ + -^-g^p^p^ + J^-g^pfjf* 
Ac w s w lui^^ Trig^i) ^£j(i) 



m BW m BW m BW 

m B(i> 



^QM^f 3 " 4 - (11) 

where we have pulled out a common factor with ay = a em /c^ and we sum over all charged 
fermions running in the loop. The fact that the WIMPs are annihilating nearly at rest (i.e., 
p/j. (mg(i),0)) also allows one to define the z-axis in the center-of-mass frame as the axis 
along which the final state particles are traveling. This simplifies things greatly as many of 
the dot products between polarization tensors and momentum vectors vanish. Accounting 



S 



for this and imposing conservation of momentum (2p = p 3 + Pa), we find that the only 
remaining tensor structures are B 2 , Ba, B G ,C\, C 2 and C3. 

Computation of the loop integrals in WIMP annihilations is complicated by the fact 
that two of the incoming momenta are nearly identical. In these configurations, usual 
approaches such as the Passarino-Veltman scheme 



26| for computing tensor integrals break 
down and one must rely on alternative schemes in order to safely compute one-loop scattering 
amplitudes. We have chosen to use the scheme introduced in Ref. [2^]. This scheme was 
most recently used to compute one-loop WIMP annihilation scattering amplitudes and we 



refer the interested reader to in Ref. 



171 ] for details. 



Squaring the amplitude and summing/averaging over all polarizations, we find that the 
cross section for + B^ — >■ 7 + Z takes the form: 

2 



(<T 7 ZV) 



1 2 2 

9 I671-4 s 2 w 



PI 



32m 2 B(1) 



Qt-AzA,k 



where (3z 



Am 



(12) 



B. Annihilation to 7 + H 



The amplitude for this process takes the form: 



M = et 1 (p 1 )e^(p 2 )e^*(p 3 )M m , 3 (p 1 ,p 2 ,P3,Pi) 



(13) 



where we have assigned the Lorentz indices and momenta as 

Bftfa) + B$(p 2 ) -+ 7ms(P3) + H( P a). 



(14) 



The amplitudes for this process involve one less vector boson, and thus have a simpler tensor 
structure than the 7Z case. In fact, due to the scalar coupling of the Higgs boson, the only 
surviving terms from the trace over the internal fermion line are proportional to Levi-Civita 
tensors e^ uXa contracted with at least one of the external momenta. Instead of expanding 
the scattering amplitude in terms of all possible permutations, we only list the surviving 
terms once the non-relativistic nature of the WIMPs and conservation of momentum are 
applied. The resulting expression takes the form: 



■H*i««»(Pl>P2,P3,P4) 



2m w s w 



Pi 

m3 B m 



Do 



m BW - 



-Pat 



+ 



D, 



-PA,ae 



HA1 



(15) 



We note that the scattering amplitude is directly proportional to the zero mode fermion 
mass rri£ (via the Yukawa coupling). As a result, only loops containing top quarks and their 
KK partners will make a significant contribution to the process. Summing and averaging 
over all polarizations, we find the cross section can be written as: 



P 2 



II 



1 2 2 

9 4nm 2 w sl \S2m 2 B(1) 



J2Q^i {yld- Y ls)^HA, 



(16) 



where j3n 



1 - 



4m 2 



(l) 



In Fig.[3j we plot the cross sections for the 77, 7Z and 7 if final states as a function of the 
WIMP mass 777,^(1) where we have assumed a universal KK fermion mass of vn^ — 1.1 x Tn B (v>. 
From this figure, we see that the 7Z cross section is roughly 10% of the 77 cross section 
and is in agreement with the estimation of Ref. [l^| while the 7 if cross section is even more 
suppressed. 



V. GAMMA RAY FLUX AND ENERGY SPECTRUM 

The differential flux of photons predicted to be observed at an angle ip from the direction 
to the galactic center (GC) is 



d$ 7 
dVtdE 



r Q p 2 Q (IN., 

:7im 2 B{1) dE 



with 



diV 7 

He 



r ds 


~p[r(s,ip))~ 


/l.O.S. r 


. P& . 






)f dE ' 





(17) 
(18) 



where the index / denotes the annihilation channels with one or more photons in the final 
state, (av) j is the corresponding cross-section and dNf /dE is the (normalized) photon spec- 
trum per annihilation. Furthermore, p(x), p & = 0.3 GeV/cm 3 and r = 8.5 kpc respectively 
denote the dark matter density at a generic location x with respect to the GC, its value at 
the solar system location and the distance of the Sun from the GC. Finally, the coordinate 
s runs along the line of sight, which in turn makes an angle ip with respect to the direction 
of the GC. 
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FIG. 3: The annihilation cross sections for the 77, jZ and final states. 



Equation ffTTj) allows one to separate the factors contributing to the predicted incoming 
gamma ray flux. In particular, while the dN^/dE term is sensitive only to the particle 
physics under consideration, the remaining factors are sensitive to the modeling of the halo 
density profile p(x). For a given dark matter model, these latter factors are the main source 
of uncertainty in the prediction of the detectability of a dark matter signal. To proceed 
further, we define with J the dimensionless integral along the line-of- sight appearing in 
Eq. (I17p and with J its average value computed for a solid angle AQ centered on the GC 

2 



J 

J (Ml) 



r ds 




)]] 


/l.O.S. r 


. Pq 





dn. 



(19) 



M JAQ 

The J factor thus defines the normalization of the gamma ray flux signal and allows one to 
quantify the impact of astrophysical uncertainties due to the lack of knowledge of the halo 
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density profile. In what follows we will consider two models for p: the Navarro Frenk and 
White (NFW) profile and the "adiabatically contracted" profile. 

The Navarro Frenk and White (NFW) profile has been shown to fit reasonably well the 



results from recent high-resolution numerical simu 
benchmark for indirect dark matter searches 



ations. This density profile is often used 



Pnfw(^) 



Ps 



2 ■ 



(20) 



Modifications of the above pro file on very small scales have been observed in the most 



recent simulations. While Ref. 
approximated with r~ 1,2 cusps, Ref. 



29| shows that the innermost regions of DM halos are better 



30] points out that the analytic form that provides an 



optimal fit to the simulated halos is the so-called "Einasto profile" 31] 



p(r) = p exp 



2 

a 



(21) 



which is shallower than NFW at very small radii 2 . 

The above results have been derived in the framework of simulations containing only 
dark matter particles, which interact gravitationally. Baryons however are expected to have 
a non-negligible effect on the small scale dynamics of galaxies. In particular, due to the 
dissipative nature of the baryonic fluid, when baryons lose energy and contract, this affects 
the gravitational potential experienced by the dark matter. In the "adiabatic compression" 



scenario 



32] , the baryons contraction is quasi-stationary and spherically symmetric. Startin 



from an initial NFW profile, the final slope in the innermost regions becomes r -15 331-136]. 

Assuming AQ = 10~ 5 sr. (corresponding to the angular resolution of the HESS and 
Fermi LAT 7-ray experiments), in Tab. [I] we show the value of J obtained for two profiles: 
the benchmark NFW profile and the "adiabatically contracted" profile, with the same pa- 



rameters as the profile labelled 'A' in Ref. 



36] 3 . Table [I] explicitly shows the extent to which 



the present uncertainty in the knowledge of the halo density profile turns into uncertainty 
in the predicted gamma ray flux from the GC. 



2 The values assumed for the parameters are in this case a = 0.17 and R = 20 kpc. 

3 Note tha t our values of J for the NFW and "Adiabatic" profiles are slightly different from the values in 
Ref. 



361 ] , due to the fact that the NFW profile was there approximated as a simple r 1 power- law from 



the galactic center to the location of the Sun. 
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Model 


J(1(T 5 ) 


NFW 


1.5 x 10 4 


Adiabatic 


4.7 x 10 7 



TABLE I: Value of J(10 5 ) for two dark matter density profiles. 



Recently, a large suite of cosmological simulations aimed to study the assembly of Milky 
Way like galaxies were carried out using the Adaptive Mesh Refinement (AMR) code RAMSES, 
which includes treatment of dark matter, gas and stars 37| . The modeling included real- 
istic recipes for star formation, supernova feedback (SNII and SNIa), stellar mass loss, gas 
cooling/ heating and metal enrichment. It was found that adopting different prescriptions 
for the star formation history results in a strong difference not only in galactic disk size and 
concentration, but also on DM contraction, therefore dramatically affecting the DM profile 
in the innermost regions of galactic halos. 

These simulations were shown to bracket the so-called Kennicutt-Schmidt relation (that 
parametrizes the star formation rate as a function of the gas density) from the THINGS 



survey 



38 



and they are therefore representative of Milky Way-like galaxies in the local 
Universe. A comparison of the profiles obtained in these simulations with the standard 
DM-only case can be found e.g. in Ref. 39[, where it is also shown that the final profile is 
substantially different with respect to a naive implementation of adiabatic contraction. We 
also stress that DM halos in numerical simulations deviate rather strongly from spherical 
symmetry, a circumstance that introduces additional uncertainties on the estimate of the 



integral along the line of sight J 



39). 



Other processes due to the gravitational interaction of baryons (such as the "heating" of 
the DM fluid) can possibly have the opposite effect on the final DM distribution and soften 
the central cusps. Furthermore, the presence of a super-massive black hole (BH) at the GC 
inevitably affects the DM profile. The growth of the BH from an initial small seed would 



40] . While dynamical effects and DM annihilations will 



36 



41 



42( | . in some cases significant overdensities 



initially lead to a large DM "spike 
subsequently tend to destroy the spike 
can survive over a Hubble time. 

In what follows, the predicted gamma ray flux from the galactic center is obtained as- 
suming the NFW halo profile. To obtain predictions for other profiles, it is sufficient to 
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FIG. 4: The gamma ray flux as a function of the photon's energy for a WIMP of mass 1 TeV. 
Shown are three different experimental energy resolutions. 

rescale the flux by the appropriate ratio of J (which is 3.3 • 10 3 for a naive implementation 
of adiabatic contraction). 

Finally, we account for the finite resolution of the detector by convolving the unfiltered 
signal S(E) with a gaussian kernel G(E, E ), 

(22) 

where a is related to the detector's relative energy resolution £ by a = £/2.3. The signal 
Sm(E ) measured by the detector at energy E is then simply given by 

S M (E ) = J dEG(E, E ) S(E). (23) 

Combining the expectations for J with the particle physics rates and the detector re- 
sponse, we arrive at predictions for the flux and spectrum of gamma rays. In Figs. HI |5]and 
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FIG. 5: The gamma ray flux as a function of the photon's energy for a WIMP of mass 500 GeV. 
Shown are three different experimental energy resolutions. 



3 mass. 



wc 



[6j we show predictions for the flux for an NFW halo profile and three choices of LK 
Although the thermal relic density favors LKPs in the range of 500 - 900 GeV 
consider a wider range of LKP masses motivated by the possibility that the LKP is not a 
canonical thermal relic, which could allow for a wider range of masses consistent with the 
observed quantity of dark matter in the Universe. For example, if the modulus which con- 
trols the size of the extra dimension is cosmologically active, it could cause large deviations 



in the thermodynamics of the Universe from its extrapolated history (43] . Alternately, early 
production of KK gravitons which eventually decay into LKPs can serve as a nonthermal 
production mechanism for KK dark matter in an otherwise standard cosmological history 



44] 



Figures H]- [6] contain curves for three different choices of energy resolutions, ranging from 
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FIG. 6: The gamma ray flux as a function of the photon's energy for a WIMP of mass 300 GeV. 
Shown are three different experimental energy resolutions. 

the ~ 10% resolution typical of current experiments, to an aggressive 0.5% resolution which 
might be possible in future experiments. We find that at 10% energy resolution, lines in 
the 5d UED model are very difficult to distinguish from the continuum. At a 5% energy 
resolution, broad lines may appear for LKP masses around 300 GeV, slightly above the lower 
bound from colliders. At 0.5%, well separated lines for 77, jZ, and 'jH are visible for light 
LKPs, and some structure related to the jH line is visible for an LKP mass of around 500 
GeV. 

In principle, we should compare our predicted flux with gamma ray observations, since 
data are available from a variety of gamma-ray telescopes, such as the Fermi LAT and Air 
Cherenkov Telescopes like HESS and MAGIC. The comparison is however made complicated 
by the aforementioned uncertainties on the normalization of the predicted flux on one side, 
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and on the other by the difficult interpretation of astrophysical data. 

It is for instance unclear at this stage whether known astrophysical sources can account for 
the energy spectrum obtained by Fermi. A recent analysis of a 1° x 1° found that the diffuse 
gamma-ray background and discrete sources, as estimated with standard astrophysical tools, 
can account for practically all the observed flux, but there is a small excess at energies 0(1) 
GeV, which may or may not be explained with an improved model of the Galactic diffuse 
emission and with a better understanding of the systematic errors {45]. At the same time, 
there is an ongoing debate on the possible origin of the point source observed by HESS right 



at the Galactic center 



46). 



In principle, once these issues are settled, and the astrophysical backgrounds are well- 
known, one can search for a DM signal even in the case where it provides a subdominant 
contribution to the total flux. More subtle features such as subdominant gamma-ray lines 
could therefore in principle emerge in the process of data analysis, in which case one could 
try to perform the program describe above, and search for additional lines that might shed 
light on the structure of the dark sector. 



VI. CONCLUSIONS 



Extra dimensions are a fascinating possibility in the spectrum of physics beyond the 
Standard Model, one which leads to a unique mechanism to provide a stable WIMP to play 
the role of dark matter. We have examined the gamma ray line signals for UED dark matter 

n 

WIMPs, extending the results for 77 in the literature [14J to also include lines from 7Z and 
^yH final states. Our conclusions are that such lines are very challenging for the current 
generation of gamma ray observatories to resolve, but may be visible at the next generation 
of such experiments. As such experiments are designed, we hope our work will be of some 
value in evaluating their exciting physics capabilities. 
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